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Dynamical counterflow effects on vortex evolution under the truncated Gross-Pitaevskii equation 
are investigated. Standard longitudinal mutual friction effects are produced and a dilatation of 
vortex rings is obtained at large counterflow. A strong temperature-dependent anomalous slowdown 
of vortex rings is observed and attributed to the presence of thermally exited Kelvin waves. This 
generic effect of finite-temperature superfluids is estimated using energy equipartition and orders 
of magnitude are given for weakly interacting Bose-Einstein condensates and superfluid ^He. The 
relevance of thermally excited Kelvin waves is discussed in the context of quantum turbulence. 
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Quantum vortices present in superfluids interact with 
the normal fluid producing mutual friction effects that 
must be phenomenologically introduced into Landau's 
two-fluid model [ll [2]. For superfluid ^He, there is 
no generally-accepted theory of mutual friction that 
is valid over the entire temperature range [3 . For 
Bose-Einstein Condensates (BEG), the Gross-Pitaevskii 
equation (GPE) is a dynamical description that was 
thought to be valid only in the low-temperature limit 
[4]. Davis et al. [5 suggested that, when a truncation of 
Fourier modes is performed, the resulting truncated GPE 
(TGPE) can also describe the (classical) thermodynamic 
equilibrium of homogeneous BEG [5]. The TGPE was 
found to relax toward (microcanonical) equilibrium and 
a condensation transition was obtained [3 0- Vortex 
dynamics was studied within the TGPE by Berloff and 
Youd ^ who observed a dissipative contraction of vortex 
rings. 

The purpose of this Letter is to investigate mutual fric- 
tion and counterflow effects in the context of the TGPE. 
We present a stochastic algorithm that allows to effi- 
ciently generate grand canonical equilibrium states with 
non-zero momentum at given (target) values of tempera- 
ture chemical potential and counterflow. These states are 
then combined with lattices of straight vortices and vor- 
tex rings and their TGPE evolutions are monitored. Our 
main result is that, beside the phenomenologically ex- 
pected counterflow effects, the TGPE also induces a (phe- 
nomenologically) unexpected slowdown of vortex rings 
that is caused by thermally excited Kelvin waves and 
should be considered in quantum turbulence. 

The TGPE describing a homogeneous BEG of volume 
V is obtained from the GPE by truncating the Fourier 
transform of the wavefunction ?/^: = for |k| > /cmax 
[H 1^ . Introducing the Galerkin projector Vq that reads 
in Fourier space T^ciV^k] = ^(^max — |k|)V^k with 0{-) the 
Heaviside function, the TGPE explicitly reads 
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where is the number of condensed particles per unit 
volume, m is their mass and g = 47rafi^/m, with a the 
s-wave scattering length. The superfluid velocity reads 
Vg = {h/m)\/(j)^ where (j) is the phase of the (complex) 
i/j and h/m is the Onsager-Feynman quantum of velocity 
circulation around vortex lines i/j = [4 . When Eq. ([T]) 
is linearized around a constant t/j = i/jq^ the sound veloc- 
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ity is given by c = (^IV^oP/^) with dispersive effects 
taking place for length scales smaller than the coherence 
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length ^ = (fi^/2m|?/;oP^) that also corresponds to the 
vortex core size. 

Equation ([T]) exactly conserves the energy H = 

Jd^x (l^lVV^p + f [T^gIV^P]^) and the number of par- 
ticles N = J (Px\ip\'^. The momentum P = 
Y J d^x (V^VV^ — t/jS/t/j) is also conserved when standard 
Fourier pseudo- spectral methods are used, provided that 
they are dealiased using the 2/3-rule (/cmax = 2/3 x M/2 
[8^ at resolution M) [T5] . 

Microcanonical equilibrium states are known to re- 
sult from long-time integration of TGPE [SHZ]- Grand 
canonical states are given by the probability distribution 
Pst[V^] = exp[-p{H - /iAT - Vn • P)]. They allow to 
directly control the temperature (instead of the energy in 
a microcanonical framework). These states can be effi- 
ciently obtained by constructing a stochastic process that 
converges to a realization with the probability Pgt [V^] [Sj • 
This process is defined by a Langevin equation consisting 
in a stochastic Ginbzurg-Landau equation (SGLE): 
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where the white noise C(x?^) satisfies (C(x, t)C*(x', t')) = 
5(t — t')5{K — x'), 13 is the inverse temperature, ji the 
chemical potential and Vn the normal velocity. The term 
ifiVn • VV^ induces an asymmetry in the repartition of 
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sound waves and generates non-zero momentum states. 
These states do not generally correspond to a condensate 
moving at velocity Vg = Vn because Vg is the gradient 
of a phase and takes discrete values for finite size sys- 
tems. Equilibrium states states with nonzero values of 
the counterflow w = Vn — Vg are generated in this way. 

Using this algorithm in |9] the microcanonical and 
grand canonical ensembles were shown to be equivalent 
and the condensation transition reported in [5l |6] identi- 
fied with the standard second order A-transition. All the 
SGLE equilibrium used in this letter have a condensate 
at rest (vg = 0) and therefore Vn = w. 

At low-temperature the partition function Z can be 
exactly computed by the steepest-descent method [9 . In 
particular, setting Vn = (0, 0, the momentum and the 
number of particles of the equilibrium state read = 
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where 



/ [w^^l and N 

is the total number of modes and f[z] 
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(y/z). These relations furnish an explicit 
expression for the normal density pn = y 

The direct control of the counterfiow in the SGLE 
algorithm allows to obtain the temperature dependence 
of pn in the TOPE context. Low-temperature exact re- 
sults are in good agreement with SGLE data, see Figjl] 

In all the numerical simulations presented in this Let- 
ter jii is adjusted in order to fix the density p = mN/V 
to 1 and the physical constants in Eqs.Q and ([2| are 
determined by the relations ^/cmax = 1.48 and c = 2. 
The inverse temperature is normalized as /3 = M jVT 
and V = (27r)^. With these choice of parametrization 
the A-transition temperature is independent of M and its 
value is fixed to T\ = 2.48; the quantum of circulation 
him has the value 




FIG. 1: Temperature dependence of the normal density pn 
(see text). Inset: Pz as a function of at fixed temperature 
T — O.OSTa. Points: SGLE ([2| equilibration at resolution 
64^; solid lines: low-temperature exact results. 



We now turn to counterfiow effects. To wit, we use 
an array of alternate-sign straight vortices V^iattice (see 
[To]). This exact stationary solution of the GPE is ob- 
tained by a Newton method. The vortices are separated 
by a distance tt and can be considered isolated when 
^ ^ 0, as the resolution is increased. An equilibrium 
state ?/^eq is prepared using the SGLE ([2| with counter- 
fiow perpendicular to the vortices. The initial condi- 



tion = '^/^lattice X V^eq IS then cvolvcd with the TGPE. 
Figure [2ja displays 3D visualizations of the density at 
t = and t = 100 where the displacement of the lat- 
tice is apparent. The temporal evolution of the (parallel 
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FIG. 2: a-b) Density at t = 0, 100 of the lattice configuration 
(red) with T — 0.4 Ta and — 0.4. Blue clouds correspond to 
density fluctuations, c) Positions {R\\^R±_) of a single vortex 
for T = 0.2Ta, T = 0.4 Ta and Vr, = .4. Resolution 64^ 



and perpendicular to Vn) position of a vortex {R\\^R^) 
are presented on Fig|2lc for T = 0.2 Ta, T = 0.4 Ta and 
= 0.4. The counterfiow-induced vortex velocity clearly 
depends on the temperature. A perpendicular motion is 
also induced at short times. This motion has two phases: 
first an adaptation, making the lattice slightly imper- 
fect, followed by a much slower perpendicular motion. 
Observe that the imperfection of the lattice at final con- 
figurations is almost equal for the two temperatures pre- 
sented in Fig|2]c, but the parallel velocities are consider- 
ably different. The self-induced parallel velocity caused 
by the slight lattice imperfection is thus very small and 
not driving the longitudinal motion. 

We now concentrate on the measurement of R\\ for 
which the present configuration is best suited. R\\ has a 
linear behavior, that allows to directly measure the par- 
allel velocity v\\ . The temperature dependence of v\\/v^ is 
presented on Figj4]for different values of and ^. This 
behavior is consistent with the standard phenomenolog- 
ical model for the vortex line velocity vl [3 : 

vl = Vsi + as' X (vn - Vsi) - a's' X [s' X (Vn - Vsl)], (3) 

where s' is the tangent of the vortex line, Vgi = Vg + ui 
is the local superfiuid velocity with ui the self-induced 
vortex velocity and Vn the normal velocity. The mutual 
friction coefficients in Eq.Q are typically written as a = 
5pn/2p, a' = B' p^/2p where B and B' are order-one and 
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weakly temperature-dependent. Equation (|3| applied to 
a straight vortex with perpendicular to the vortex and 
Vs = yields = v\\/vn- The value of = B' pn/2p with 
B' = 0.83 is displayed on Fig|4] (bottom dashed line) and 
is in good agreement with the lattice data. 

We now turn to the interaction of vortex rings and 
counterflow. The Biot-Savart self-induced velocity of a 
perfectly circular vortex ring of radius R is given by 



h C{R/0 
2m R 



C{z) = In (8z) 
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where a is a core model-depending constant [3 . We have 
checked, using an initial data V^ring prepared by a New- 
ton method that the GPE (large R/^) ring translational 
velocity is well reproduced by Q with a = 0.615. 

Equation ([3| with perpendicular to the ring and 
Vs = yields the radial velocity R = —a{ui — v^). The 
case without counterflow {v^ = 0) was studied by Berloff 
and Youd [7] and a contraction of vortex rings compatible 
with Q was reported. To study the influence of counter- 
flow we prepare an initial condition t/j = V^ring x '^eq in the 
same way as above for the vortex lattice. The temporal 
evolution of the (squared) vortex length of a ring of initial 
radius R = 15^ at temperature T = 0.4 Ta and = 0, 
0.2 and 0.4 is displayed on FigjSja. The Berloff- Youd con- 
traction [7 is apparent in absence of counterflow (bottom 
curve). The temperature dependence of the contraction, 
related to the a coefficient in Eq.(|3|, also quantitatively 
agrees with their published results (data not shown). 




FIG. 3: a) Temporal evolution of the (squared) length of a 
vortex ring at different values of counterflow Vn (temperature 
T — 0.4 Ta and initial radius R — 15^). b-c) 3D visualization 
of vortex ring {R = 20^) and density fluctuations at t = 18, 
19, with T — 0.4 Tx and resolution 64^. Same colorbar as in 
Fig [2] Thermally-excited Kelvin waves are apparent. 



A dilatation of vortex rings is obtained (top curve on 
Figjsja). when the counterflow v^^ is large enough. Such a 
dilatation -a hallmark of counterflow effects- is expected 



|3] to correspond to a change of sign of Vn — Vgi in Eq. (|3|. 
However, the predictions of Eq.Q unexpectedly turn out 
to be quantitatively wrong. Indeed, using Eq.Q in the 
conditions of Fig|3]a one finds Vgi = ui = 0.39 which is 
significantly larger than normal velocity = 0.2 around 
which dilatation starts to take place (see middle curve 
on Figjsja). Equation ([3| prediction for the longitudinal 
velocity vl = {l—a')u\^a'vn is also unexpectedly wrong. 
Using the value of a' determined above on the vortex 
array, one finds vl 0.98iii and from Eq.Q one finds 
for Vl the value 0.38 that is larger than the measured 
value Vl = 0.23. 

This anomaly of the ring velocity vl is also present in 
the absence of counterflow ('^n = 0) where Eq.(|3]) predicts 
that a' should be given by Avl/u^ = {u^ — vl)/u{. The 
temperature dependence of /S.vl/u\ is displayed on Figj4] 
(top curve). Observe that Avl/ui is one order of mag- 
nitude above the transverse mutual friction coefficient a' 
measured on the lattice. 
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FIG. 4: Temperature dependence of counterflow-induced lat- 
tice velocity v\\/vn (bottom) and ring slowdown Avl/ui (top) 
obtained with Vn — 0. Dashed line: prediction of Eq.(|3| with 
a — 0.83/?n/2p; solid line: prediction of anomalous slowdown 
by Eq.(|6l with R — 20^ at various resolutions. 



We now relate the thermally-induced anomaly to the 
velocity induced on a vortex ring by a single a Kelvin 
wave of (small) amplitude A and (large) wavenumber 
Nk/^ttR obtained in the LIA fTTl and Biot-Savart [T2] 
frameworks. The velocity Vg, reads (see Eq. (26) of [TT]) 



(5) 



where is the (undisturbed) ring velocity Q. 

The TGPE model naturally include thermal fluctua- 
tions that excite Kelvin waves as apparent on Figj3]b-c. 
We assume that the slowing down effect of each individ- 
ual Kelvin wave is additive and that the waves populate 
all the possible modes. Kelvin waves being bending oscil- 
lations of the the quantized vortex lines their wavenum- 
ber must satisfy k < = 27r/^. The total number of 
thermally excited Kelvin waves is thus A/xeivin ^ Rk^. 

The amplitude term A^N'^/R^ in ([5| can be obtained 
by simple equipartition arguments. The energy of a 
(perfect) ring E = ^^^^^R[C{R/^) - 1], with 
the superfiuid density [3 . A Kelvin wave produces a 
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variation of the ring length AL = ttA^N^/ R. Its en- 
ergy can thus be estimated as AE = Assuming 
AE = /3~^ yields, at low temperature where ps ^ 
A^Nl/R^ = m^p-^/7T^ph'^RC{R/0 Replacing 
jR? in Eq.(|5]), the dominant effect is obtained by sum- 
ming up to A/xeivin and it finally reads: 

Avl Uj - _ p-^rri^ 
u, ~ u, ^ 7r^ph^C{R/0^' 

The thermally-induced anomalous slowdown ([6| is in 
good agreement with the TGPE data displayed on Figj4] 

We now extend (|6| in order to take into account quan- 
tum effects and estimate orders of magnitude in the phys- 
ical case of BEG and superfluid ^He. The dispersion 
relation of Kelvin waves uj{k) = ^k'^C{R/S,) [11 im- 
plies (using the relation Jiw{keq) = = k^T) that 
Kelvin waves are not in equipartition for wavenumbers 
k> kec^ = {2mkBT/h^C{R/i)Y^'^, as (like in blackbody 
radiation) quantum effects are relevant in this range. 

For weakly-interacting BEG with mean inter-atomic 
particle distance I ~ IV^ol"^^^ satisfying a <C ^ <C ^ 
the condensation temperature is Tx ~ /k^ml'^. For 
T > T*, where T"" /Tx - C{R/i) fji^ <Cl, it is straight- 
forward to show that /Cgq > % and therefore that (|6| di- 
rectly applies and reads Av^/^i ^ {l/i){T /TxC{R/i)). 
For T < T* ^ k^ must be replaced by /Ceq in formula (|6| 
and the slowdown becomes Avl/u; ^ (T /Tx C{R/£/)Y'^. 

At zero-temperature it is natural to suggest that the 
quantum fluctuations of the amplitudes of Kelvin waves 
produce an additional effect. This effect can be estimated 
by using AE = fux){k)/2. It is radius-independent and 
of order Avl/ui ^ (^0^ (^^^ 0)- interesting to 

note that a drift effect of order {i/^)^ is also obtained 
by balancing with a Magnus force [3] the Roberts and 
Pomeau [13] Gasimir-like force due to the scattering of 
(density) zero-temperature quantum fluctuations. 

In a low-T physical BEG, with quantum distribution of 
sound waves, pn/ p ^ {T/Tx)^ PJ and the standard effects 
(|3| are of order {T/TxY. Thus the new effect should 
dominate in this limit. In the case of superfluid ^He 
the GPE description is only expected to give qualitative 
predictions Nevertheless the new effect should also 
be dominant at low- temperature. 

Thermally excited Kelvin waves can also be relevant in 
the broader context of quantum turbulence where Kelvin 
waves are excited at low temperature by vortex reconnec- 
tion and their energy finally decays into sound waves [14]. 
The phenomenological speed of vortex rings (Eqs.([3| and 
Q) is routinely used in this context, clearly an overesti- 



mation as these rings are thermally perturbed and hence, 
in reality, slower. In a dilute gas of vortex rings this slow- 
down will increase the time between collisions and inhibit 
reconnection by a factor given by ([6|. In the case of a 
dense vortex tangle their effect is more difficult to predict 
but should be, in the future, studied using the TGPE. 

In summary we obtained and measured standard coun- 
terflow mutual friction effects within the TGPE. Our 
main result is that vortex rings are decelerated by ther- 
mal fluctuations of Kelvin waves and that these fluc- 
tuations, generic of finite-temperature superfluids, pro- 
duce an experimentally-testable effect that dominates the 
standard effects at low-temperature. 

We acknowledge useful scientific discussions with G. 
Barenghi. The computations were carried out at IDRIS 
(GNRS) and visualizations used VAPOR [l^. 
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